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Abstract 



We investigate the possibility of detecting the Higgs bosons predicted in the 
Minimal Super symmetric extension of the Standard Model (h° , H° , A , H^) , 
with the reactions e + e~ — > bbh (H°,A ), and e + e~ — ► t~ v T H + ,t + v t H~ , 
using the helicity formalism. We analyze the region of parameter space (m^o — 
tan/3) where h°,H°,A° and could be detected in the limit when tan/3 
is large. The numerical computation is done considering two stages of a 
possible Next Linear e + e~ Collider: the first with y/s = 500 GeV and design 
luminosity 50 fb^ 1 , and the second with yfs = 1 TeV and luminosity 100-200 
fb- 1 . 

PACS: 14.80.Cp, 12.60.Jv 
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I. INTRODUCTION 



Higgs bosons M play an imporant role in the Standard Model (SM) 0; they are respon- 
sible for generating the masses of all the elementary particles (leptons, quarks, and gauge 
bosons). However, the Higgs-boson sector is the least tested one in the SM. If Higgs bosons 
are responsible for breaking the symmetry from SU(2)l x U(1)y to U(1)em, it is natural 
to expect that other Higgs bosons are also involved in breaking other symmetries at the 
grand-unification scale. One of the more attractive extensions of the SM is Supersymmetry 
(SUSY) P|, mainly because of its capacity to solve the naturalness and hierarchy problems 
while maintaining the Higgs bosons elementary. 

The minimal supersymmetric extension of the Standard Model (MSSM) doubles the 
spectrum of particles of the SM and the new free parameters obey simple relations. The 
scalar sector of the MSSM |4| requires two Higgs doublets, thus the remaining scalar spec- 
trum contains the following physical states: two CP-even Higgs scalar (h° and if ) with 
m h o < m H o, one CP-odd Higgs scalar (A ) and a charged Higgs pair (if 1 * 1 ), whose detection 
would be a clear signal of new physics. The Higgs sector is specified at tree-level by fixing 
two parameters, which can be chosen as the mass of the pseudoscalar m^o and the ratio of 
vacuum expectation values of the two doublets tan/3 = V2/V1, then the mass m h o, m H o and 
rriH± and the mixing angle of the neutral Higgs sector a can be fixed. However, since radia- 
tive corrections produce substantial effects on the predictions of the model 0, it is necessary 
to specify also the squark masses, which are assumed to be degenerated. In this paper, we 
focus on the phenomenology of the neutral CP-even and CP-odd scalar (h°, H°, A ) and 
charged {H^. 

The search for these scalars has begun at LEP, and current low energy bound on their 
masses gives m h o, uiao > 90 GeV and ra#± > 120 GeV for tan (3 > 1 1|. At e + e~ colliders the 
signals for Higgs bosons are relatively clean and the opportunities for discovery and detailed 
study will be excellent. The most important processes for the production and detection of 
the neutral and charged Higgs bosons h°, H°, A and H^, are: e + e~ — > Z* — ■> h°, H° + Z°, 
e+e" -> Z* -> h°, H° + A°, e + e~ -> uu + W + *W~* -> uv+h°, H° (the later is conventionally 
referred to as WW fusions), and e + e~ — > H + H~ JJJ; precise cross-section formulas appear 
in Ref. ||. The main decay modes of the neutral Higgs particles are in general bb decays 
(~ 90%) and t + t~ decays (~ 10%) which are easy to detect experimentally at e + e~ colliders 
IPHlljl • Charged Higgs particles decay predominantly into rv T and tb pairs. 

The Z°h° production cross-section contains an overall factor sin 2 (/? — a) which suppress 
it in certain parameter regions (with m^o < 100 GeV and tan (3 large); fortunately the A°h° 
production cross-section contains the complementary factor cos 2 (f3 — a). Hence the Z°h° 
and A°h° channels together are well suited to cover all regions in the (m A o — tan 0) plane, 
provided that the cm. energy is high enough for Z°h° to be produced through the whole 
m^o mass range, and that an adequate event rate can be achieved. These conditions are 
already shown to be satisfied P~ P11 f° r = 500 GeV with assumed luminosity 10 fb^ 1 , 
as is expected to be the case of the Next Linear e + e~ Collider (NLC). 

In previous studies, the two-body processes e + e~ — ► h°(H°)+Z° and e + e~ — > h°(H°)+A° 
have been evaluated f| extensively. However, the inclusion of three-body process e + e~ — > 
h°(H°) + bb and e + e~ — >• A + bb ]I2] at future e + e~ colliders energies [|I3|-|T5| is necessary 
in order to know its impact on the two-body mode processes and also to search for new 
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relations that could have a cleaner signature of the Higgs bosons production. 
In the other hand the decay modes of the charged Higgs bosons determine the signatures in 
the detector. If m H ± > m t + m b , the dominant decays modes are tb, ib and t + u t , t~v t . In 
some part of the parameter space also the decay H + — > W + h° is allowed. If m H ± < m t + m b , 
the charged Higgs boson will decay mainly into t + u t , t~v t . 

For mpfi < mzo, and if 50 events criterion are adecuate, the H + H~ pair production will 
be kinematically allowed and easily observable |Q-|ri|,|T^-[TB| . For m^o > 120 GeV, e + e~ — > 
H + H~ must be employed for detection of the three heavy Higgs bosons. Assuming that 
SUSY decays are not dominant, and using the 50 event criterion H + H~ can be detected up 
to m H ± = 230 GeV @^,[lj|4l]| , assuming = 500 GeV. 

The upper limits in the H + H~ mode are almost entirely a function of the machine energy 
(assuming an appropriately higher integrated luminosity is available at a higher \^s). Two 
recent studies [|l9l|20l show that at yfs = 1 TeV, with an integrated luminosity of 200 fb" 1 , 
H + H~ detection would extended to m^o ~ m H ± ~ 450 GeV even if substantial SUSY 
decays of these heavier Higgs are present. 

In the present paper we study the production of SUSY Higgs bosons at e + e~ colliders. 
We are interested in finding regions that could allow the detection of the SUSY Higgs bosons 
for the set parameter space (m^o — tan/5). We shall discuss the neutral and charged Higgs 
bosons production bbh°(H°, A ), and t~v t H + , t + v t H~ in the energy range of a future e + e~ 
colliders fl3Hl5f for large values of the parameter tan (3, where one expects to have a high 
production. Since the coupling h°bb is proportional to sin a/ cos/5, the cross-section will 
receive a large enhancement factor when tan/3 is large. Similar situation occurs for H°, 
whose coupling with bb is proportional to cos a/ cos (3. The couplings of A with bb and of 
with t u T , t + v t are directly proportional to tan f3, thus the amplitudes will always grow 
with tan f3. We consider the complete set of Feynman diagramas at tree-level and use the 
helicity formalism |21~H 27l] for the evaluation of the amplitudes. The results obtained for 
the three-body processes are compared with the dominant two-body mode reactions for the 
plane (m^o — tan/5). Succintly, our aim in this work is to analyze how much the results 
of the Bjorken Mechanism [Fig. 1, (1.4)] would be enhanced by the contribution from the 
diagrams depicted in Figs. 1.1-1.3, 1.5 and 1.6 in which the SUSY Higgs bosons are radiated 
by a b(b) quark. For the case of the charged Higgs bosons the two-body mode [Figs. 3.1 and 
3.4] would be enhanced by the contribution from the diagrams depicted in Figs. 3.2, 3.3, 
and 3.5, in which the charged Higgs boson is radiated by a t~v t (t + v t ) lepton. 

Recently, it has been shown that for large values of tan f3 the detection of SUSY Higgs 
bosons is possible at FNAL and LHC p8| . In the papers cited in Ref. |28j the authors 
calculated the corresponding three-body diagrams for hadron collisions. They pointed out 
the importance of a large bottom Yukawa coupling at hadron colliders and showed that the 
Tevatron collider may be a good place for detecting SUSY Higgs bosons. In the case of the 
hadron colliders the three-body diagrams come from gluon fusion and this fact makes the 
contribution from these diagrams more important, due to the gluon abundance inside the 
hadrons. The advantage for the case of e + e~ colliders is that the signals of the processes are 
cleaner. 

This paper is organized as follows. We present in Sec. II the relevant details of the 
calculations. Sec. Ill contains the results for the processes e + e~ — > bbh°(H°,A°) and 
e + e~ — > t~u t H + , t + v t H~ at future e + e~ colliders. Finally, Sec. IV contains our conclusions. 
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II. HELICITY AMPLITUDE FOR HIGGS BOSONS PRODUCTION 



When the number of Feynman diagrams is increased, the calculation of the amplitude is 



a rather unpleasant task. Some algebraic forms [2j| can be used in it to avoid manual calcu- 
lation, but sometimes the lengthy printed output from the computer is overwhelming, and 
one can hardly find the required results from it. The CALKUL collaboration |$0| suggested 
the Helicity Amplitude Method (HAM) which can simplify the calculation remarkably and 
hence make the manual calculation realistic. 

In this section we describe in brief the evaluation of the amplitudes at tree-level, for 
e+e" -» bbh°(H°,A ) and e + e" r-jy T H + \t + u t H~ using the HAM [0-g7|. This method 
is a powerful technique for computing helicity amplitudes for multiparticle processes involv- 
ing massles spin-1/2 and spin-1 particles. Generalization of this method that incorporates 
massive spin-1/2 and spin-1 particles, is given in Ref. [27]. This algebra is easy to program 
and more efficient than computing the Dirac algebra. 

A Higgs boson h°, H°, A , and H ± can be produced in scattering e + e~ via the following 
processes: 

e+e~ -> bbh°, (1) 
e+e" -> bbH°, (2) 
e+e" -> bbA°, (3) 
e+e- -> r-u T H + , t + v t H-. (4) 

The diagrams of Feynman, which contribute at tree-level to the different reaction mech- 
anisms, are depicted in Figs. 1-3. Using the Feynman rules given by the Minimal Supersym- 
metric Standard Model (MSSM), as summarized in Ref. |J, we can write the amplitudes for 
these reactions [|12|]. For the evaluation of the amplitudes we have used the spinor- helicity 
technique of Xu, Zhang and Chang [^2| (denoted henceforth by XZC), which is a modifi- 
cation of the technique developed by the CALKUL collaboration f3~0fl . Following XZC, we 
introduce a very useful notation for the calculation of the processes (l)-(3) [|12|] and (4). The 



complete formulas of the processes (l)-(3) are given in Ref. [0. For the case of process (4), 
we present the relevant details of the calculations. 
Let us consider the process 

e"(pi) + e + (p 2 ) -> {T~(k 2 ) + u T (k 3 ) + H + (h), r + (k 2 ) + u T (k 3 ) + H'ih)}, (5) 

in which the helicity amplitude is denoted by Ai[X(e~), A(e + ), A(r T ), X(u T )]. 

Due to charge invariance, the cross-sections for the production processes e~e + — > r~v T H + 
and e~e + — > t + v t H~ are exactly the same. One should thus calculate the cross-section of 
one of the two reactions and then multiply by a factor of two to take into account the charge 
conjugate final state. This will enormously simplify the Feynman diagrams as well as the 
amplitudes of transition. The Feynman diagrams for this process are shown in Fig. 3. From 
this figure it follows that the amplitudes that correspond to each graph are 

Mt = -iC x P H -{k 2 + k 3 )P z ( Pl + P2 )T U 
M 2 = iC 2 P T (h + k 3 )P z (pi+p 2 )T 2 , 

M 3 = -tCMh + k 2 )P z ( Pl + p 2 )T 3 , (6) 
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where 



M 4 = -iC 4 P H -(k 2 + k 3 )P^ Pl +p 2 )T 4 , 
M 5 = iC 5 P T (h + A; 3 )P 7 (pi +p 2 )T 5 , 



g 3 m T cos26 w 

Ci = y= tanp 5-7—- 

16V2 cos^ 6V 

Q 3 Tfl T 1 

C 2 = 1= — — tan f3 



32^/2 m w cos 2 Ow' 
C 3 = C 2 , (7) 



,3 



G4 = — 7= tan p sin 

C5 = C4, 



while the propagators are 

Pz°(pi +P2) 
P H ±(k 2 + k 3 ) 



(s - m| ) + im z oT z o 
(s - m| ) 2 + (m z oT z o) 2 ' 

(2k 2 ■ k 3 - m 2 H± ) + im H T H ± 
(2k 2 -h-m 2 H± Y + (m H ±T H ±y 



where s = (pi + p 2 ) 2 and the corresponding tensors are 

Tf = u(A; 2 )(l - l5 )v{h)v{p 2 ){h -h~ h)« ~ <75)w(pi), 

T 2 M = - a* TB )(& + 3 )(1 - 75)^3)^2)7^ - <lMPi), 

T» = u(k 2 )(l - 75)^ + ^ 2 )7m« - <7 5 )^(fc3)^(p 2 )7 M « - < 75 )u(pi), (9) 

T 4 M = u(fc 2 )(l - 75)^3)^2)^1 - ^ 2 - # 3 )w(Pi), 

n = Hhh^i + &X 1 ~ lMk3H P2 )^u( Pl ). 

In fact, we rearrange the tensors T s in such a way that they become appropriate to a 
computer program. Then, following the rules from helicity calculus formalism [[H|-|27j and 
using identities of the type 

{u\{pi)l tl ux{p2)}l^ = 2ux(p 2 )ux(pi) + 2u^\(p 1 )u^\(p 2 ), (10) 

which is in fact the so called Chisholm identity, and 

i> = u x (p)u x (p) +u_x{p)u-x(p), (11) 

defined as a sum of the two projections u\{p)u\{p) and U-\{p)u-\{p). 
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The spinor products are given by 



s(Pi,Pj) = u + (pi)u^(pj) = -s(pj,pi), 
t(Pt,Pj) = u-(pi)u + (pj) = [s(pj,Pi)]*. 



(12) 



Using Eqs. (10)-(12), which are proved in Ref. [27], we can reduce many amplitudes to 
expressions involving only spinor products. 

Evaluating the tensors of Eq. (9) for each combination of (A, A ) with A, A = ±1 one 
obtains the following expressions: 



Mi(+, 


+) 


= F 1 f+' + s(k 2 ,k 3 )[s(p 2 ,k 1 )t(k 1 


,Pi) ~ 


s(p 2 , k 2 )t(k 2 


Pi) ~ 


s(P2,h)t(k 3 ,p 1 )], 




Mi{-, 


+) 


= F 1 fP + s(k 2 ,k 3 )[t(p 2 ,k 1 )s(k 1 


,Pi) - 


t(p 2 , k 2 )s(k 2 


Pi) ~ 


t(p2,k3)s(k a ,p 1 )], 


(13) 


M 2 (+. 


+) 


= F 2 f 2 ' + s(k 2 ,p 2 )t(p 1 , k l )s(k l , 


h), 










M 2 {-, 


+) 


= F 2 f2 ,+ s{k 2 ,pi)t(p 2 , k x )s{ku 


h), 








(14) 


M 3 (+, 


+) 


= F 3 f 3 ' + s(k 2 , k 1 )t(k 1 ,p 1 )s(p 2 , 


h), 










M 3 (-, 


+) 


= F 3 f-P + s(k 2 , ki)t(ki,p 2 )s(pi, 


h), 








(15) 


■M 4 (+, 


+) 


= F 4 s(k 2 , k 3 )[s(p 2 , 


- s(p 2 


, k 2 )t(k 2 ,pi) 


- s{p 2 






Ma{- 


+) 


= F 4 s(k 2 , k 3 )[t(p 2 , k 1 )s(k 1 ,pi) 


-t(p 2 


k 2 )s(k 2 ,pi) 


-t(p 2 


k 3 )s(k3,p 1 )], 


(16) 


M 5 (+. 


+) 


= F 5 s(k 2 ,p 2 )t(pi, h)s(ki, k 3 ), 












M 5 (- 


+) 


= F 5 s(k 2 ,pi)t(p 2 , h)s(ki, fc 3 ), 










(17) 



where 



F 1 = -2iC x P H ±{k 2 + k 3 )P z o( Pl +p 2 ), 
F 2 = 4iC 2 P T (h + k 3 )P z o( Pl +p 2 ), 
F 3 = -AtC 3 P u (h + k 2 )P z o( Pl +p 2 ), 
F 4 = -2iC 4 P H ±{k 2 + fc 3 )P 7 (pi +p 2 ), 
F 5 =4zC 5 P T (A; 1 + A ; 3)P 7 (pi+p 2 ), 



(18) 



and 



Here, vi 



ft + = 


«-<), 


fl' + = 


« + <), 


ft> + = 




fr + = 


my - «n 


ft + = 




fs' + = 


(v z + a z u )(v z + a z e 



•1 + 4 sin 2 9w, a z e = — 1, v* = 1 and a z v = 1, according to the experimental data 



After the evaluation of the amplitudes of the corresponding diagrams, we obtain the cross- 
sections of the analyzed processes for each point of the phase space using Eqs. (13)-(17) 
by a computer program, which makes use of the subroutine RAMBO (Random Momenta 
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Beautifully Organized). The advantages of this procedure in comparison to the traditional 



"trace technique" are discussed in Refs. piH27 



We use the Breit-Wigner propagators for the Z°, h°, H°, A and H ± bosons. The mass 
of the bottom (m b « 4.5GeV) the mass (M z o = 91.2GeV) and width (V z o = 2.4974GeV) 
of Z° have been taken as inputs; the widths of h°, H°, A and are calculated from the 
formulas given in Ref . || . In the next sections we present the numerical computation of the 
processes e + e~ — ► bbh, h = h°, H°, A and e + e~ — > t^v t H + , r + v T H~ . 



III. DETECTION OF MSSM HIGGS BOSONS AT FUTURE 
POSITRO-ELECTRON COLLIDERS ENERGIES 

In this paper, we study the detection of neutral and charged MSSM Higgs bosons at e + e~ 
colliders, including three-body mode diagrams [Figs. 1.1-1.3, 1.5, and 1.6; Figs. 2.1-2.3, 2.5 
and 2.6; Figs. 3.2, 3.3, and 3.5] besides the dominant mode diagram [Fig. 1.4; Fig. 2.4; 
Figs. 3.1, and 3.4] consider two stages of a possible Next Linear e + e~ Collider: the first 
with y/s = 500 GeV and design luminosity 50 fb~ l , and the second with y/s = 1 TeV and 
luminosity 100-200 ]b~ x . We consider the complete set of Feynman diagrams (Figs. 1-3) 
at tree-level and utilize the helicity formalism for the evaluation of their amplitudes. In the 
next subsections, we present our results for the case of the different Higgs bosons. 



A. Detection of h° 

In order to illustrate our results on the detection of the h° Higgs boson, we present graphs 
in the parameters space region (m^o — tan/5), assuming m t = 175 GeV, M~ = 500 GeV and 
tan/3 > 1 for NLC. Our results are displayed in Fig. 4, for e + e~ — > (A , Z ) + h° dominant 
mode and for the processes at three-body e + e~ — > bbh . 

The total cross-section for each contour is 0.03 pb, and 0.01 pb, which gives 1500 events, 
and 500 events to an integrated luminosity of £ = 50 fb" 1 . We can see from this figure, 
that the effect of the reaction bbh is not more important that (A , Z°) + h°, for most of 
the (m^o — tan ff) parameter space regions. Nevertheless, there are substantial portions of 
parameter space in which the discovery of the h° is not possible using either (A , Z°) + h° 
or bbh . 

For the case of y/s = 1 TeV, the results of the detection of the h are shown in Fig. 
5. It is clear from this figure that the contribution of the process e + e~ — > bbh becomes 
dominant, namely e + e~ — > (A , Z°) + h° is small in all parameter space. However, they 
could provide important information on the Higgs bosons detection. For instance, we give 
the contours for the total cross-section to, say 0.01 pb, 0.005, and 0.003 pb for both processes. 
These cross-sections give 1000 events, 500 events, and 300 events in total to a integrated 
luminosity of C = 100 fb' 1 . While for £ = 200 fb' 1 the events number is 2000, 1000, and 
600, respectively, then it will be detectable the h° at future e + e~ colliders. 
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B. Detection of H° 

To illustrate our results regarding the detection of the heavy Higgs bosons H°, we give the 
contours for the total cross-section for both processes e + e" — > (A , Z°) + H°, e + e~ — > bbH° 
in Fig. 6 for y/s = 500 GeV and C = 50 fb^ 1 . The contours for this cross-section are 0.01 
pb, 0.001 pb and 0.0001 pb for both reactions (A , Z°) + H° and M°. The number of events 
corresponding to each contour are 500, 50 and 5, respectively. 

Our estimate is that if more than 100 total events are obtained for a given process 
(A , Z°) + H° or bbH° then the Higgs boson H° can be detectable. 

For the case of y/s = 1 TeV, the results on the detection of the H° are show in Fig. 7. 
The events number for each contour is 1000, 100, and 10 for C = 100 fb' 1 and 2000, 200, 
20 for C = 200 /ft -1 . 

The effect of incorporate bbH° in the detection of the Higgs boson H° is more important 
than the case of two-body mode (A , Z°) + H°, because bbH° cover a major region in the 
parameters space {m A o — tan/3). The most important conclusion from this figure is that 
detection of all of the neutral Higgs bosons will be possible at Next Linear e + e" Collider. 

C. Detection of A° 

For the pseudoscalar A , it is interesting to consider the production mode in bbA°, since 
it can have large a cross-section due to the fact that the coupling of A with bb is directly 
proportional to tan/3, thus will always grow with it. In Fig. 8, we present the contours of 
the cross-sections for the process of our interest bbA°, at y/s = 500 GeV and C = 50 /ft -1 . 

We display the contour lines for a = 0.01,0.003,0.001, showing also the regions where 
the A can be detected. These cross-sections give a contour of production of 500, 150 and 
50 events. It is clear from this figure that very high experimental and analysis efficiencies 
are necessary for detecting the Higgs boson A . 

On the other hand, if we focus the detection of the A at y/s = 1 TeV and an integrated 
luminosity of C = 100 fb -1 , the panorama for its detection is more extensive. The Fig. 
9 shows the contours lines in the plane (m A o — tan/3), to the cross-section of bbA°. The 
contours for this cross-section correspond to 300, 100 and 10 events. While for C = 200 
fb^ 1 we have 600, 200, and 20 events respectively. It is estimated that if more than 100 
total events are obtained for bb A , then it is possible to detect the A . 

D. Detection of H± 

Our results for the H + H~ scalars are displayed in Fig. 10, 11 for e + e _ — > H + H~ 
dominant mode and for the processes at three-body e + e~ — > t~u t H + ,t + v t H~ . 

The total cross-section for this reaction with y/s = 500 GeV and C = 50 fb~ x are shown 
in Fig. 10 for each contour with 0.01, 0.001, and 0.0001 pb, which gives 500 events, 50 events, 
and 5 events, respectively. We can see from this figure that the effect of the reactions r~i' T H + 
and r + u T H~ is slightly more important than H + H~ for most of the (m^o —tan f3) parameters 
space regions. Nevertheless, there are substantial portions of parameters space in which the 
discovery of the H* 1 is not possible using either H + H~ or t~v t H + and t + v t H~. 



8 



In both cases the curves with values of 0.01 pb, 0.001 pb, and 0.0001 pb give 1000, 100, 
and 10 events to an integrated luminosity of C = 100 fb^ 1 . Meanwhile, for C = 200 fb^ 1 
we have 2000, 200, and 20 events. These cross-sections are small, however, it is precisely 
in this curve where the contribution of the processes at three-body is notable. The most 
important conclusion from this figure is that detection of the charged Higgs bosons will be 
possible at future e + e" colliders. 

IV. CONCLUSIONS 

In this paper, we have calculated the production of the neutral and charged Higgs bosons 
in association with 6-quarks and with ru T leptons via the processes e + e~ — > bbh, h = 
h°,H°,A° and e + e~ — > t~v t H + ,t + v t H~ using the helicity formalism. We find that these 
processes could help to detect the possible neutral and charged Higgs bosons at energies of 
a possible Next Linear e + e~ Collider when tan/5 is large. 

In summary, we conclude that the possibilities of detecting or excluding the neutral and 
charged Higgs bosons of the Minimal Supersymmetric Standard Model (h°, H°, A , H^) in 
the processes e + e~ — > bbh, h = h°,H°,A° and e + e _ — > t~v t H + ,t + v t H~ are important 
and in some cases are compared favorably with the dominant mode e + e~ — > (A , Z°) + h, 
h = h°,H°,A° and e + e~ — > H + H~ in the region of parameters space (m^o — tan/5) with 
tan (3 large. The detection of the Higgs boson will require the use of a future high energy 
machine like the Next Linear e + e~ Collider. 
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FIGURE CAPTIONS 



Fig. 1 Feynman Diagrams at tree-level for e + e — > bbh°. For e + e — > frfrif one has to 
make only the change sin a/ cos/5 — > cos a/ cos/5. 

Fig. 2 Feynman Diagrams at tree-level for e + e~ — > 66 A . 

Fig. 3 Feynman Diagrams at tree-level for e + e~ — > r~i? T H + ,t + u t H~ . 

Fig. 4 Total cross-sections contours in (m A o — tan /?) parameter space for e + e~ — > (A , Z°) + 
/i° and e + e~ — > with y/s = 500 GeV and an integrated luminosity of C = 50 /6 _1 . We 
have taken m t = 175 Gel 7 and M~ = 500 GeV and neglected squark mixing. 

Fig. 5 Total cross-sections contours for y/s = 1 TeV and C = 100, 200 fb' 1 . We have 
taken m t = 175 GeV, M~ = 500 GeV and neglected squark mixing. We display contours 

for e + e~ — > (A , Z°) + /i° and e + e _ — > M>/i°, in the parameters space (m A o — tan/3). 
Fig. 6 Same as in Fig. 4, but for e+e~ -> (A , Z°) + H° and e+e" -> fefcif . 
Fig. 7 Same as in Fig. 5, but for e+e" -> (A , Z°) + #° and e+e" -> &Mf °. 
Fig. 8 Same as in Fig. 4, but for e + e~ — > bbA°. 
Fig. 9 Same as in Fig. 5, but for e + e _ — > &L4°. 

Fig. 10 Same as in Fig. 4, but for e + e~ — > H + H and e + e _ — > t~v t H + ,t + v t H~ . 
Fig. 11 Same as in Fig. 5, but for e + e~ — > H + H and e + e~ — > t~v t H + ,t + v t H~. 
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